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Abstract 

Dynamical systems of the gauge glass are investigated by the method of the gauge transforma- 
tion. Both stochastic and deterministic dynamics are treated. Several exact relations are derived 
among dynamical quantities such as equilibrium and nonequilibrium auto-correlation functions, 
relaxation functions of order parameter and internal energy. They provide physical properties in 
terms of dynamics in the SG phase, a possible mixed phase and the Griffiths phase, the multicrit- 
ical dynamics and the aging phenomenon. We also have a plausible argument for the absence of 
re-entrant transition in two or higher dimensions. 
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I. INTRODUCTION 



The physics of disordered systems has been one of most fascinating subjects for theorists 
and experimentalists. The gauge glass is a typical example in this category and has attracted 
much attentions. The model describes thermodynamics of various systems such as disordered 
magnets with random Dzyaloshinskii-Moriya interaction [lj], Josephson-j unction arrays with 
positional disorder in a magnetic field 3|, vortex glasses 0], and crystal systems on disordered 
substrates 4J. In three dimensions, the spin glass (SG) transition 74] for strongly disordered 
regime has been confirmed in the gauge glass systems by numerical simulations 
y,13, 11, 12! k| and renormalization group (RG) analyses jl^ 
experimental observations ^6, 17, 



which is consistent with 



In two dimensions, there is a controversy about the existence of glass-like phase in strongly 
disordered regime. Although the long-range SG order is denied in two dimensions Q], there 
is a possibility of a quasi-long-range order in which the SG correlation decays according 



(KT) phase |22j, 



to a power law 2Q , |2J| like the ferromagnetic (FM) correlation in the Koster 



231 ] . Some numerical simulations supported the latter case 



while experimental observation [28| and numerical simulations |9j, [10 



13 



29 



24 



itz-Thouless 
26, E|, 



25 



30( deny it. In 



weakly disordered regime, the instability of the KT p 
by the perturbation expansion and the RG analysis 



ras e ag ainst small disorder is suggested 



simulations 



31. 



32L while it is denied by numerical 
3 1 1, I35I 36 , [37 1 . There is another 



and other RG analyses 
controversy in two dimensions about the existence of re-entrant transition from the KT 
)hase to the non-KT one. Earlier works with real space RG analyses suggest re-entrance 
ll, l2l |38L |39j , however the analysis has been modified and provides the absence of it 

-a 



36, 37]. The same results are obtained by Monte Carlo simulations |4J,|9 



11 



25, 




the RG analyses |3j 

Since randomness and frustration make it difficult to examine short-range systems analyt- 
ically as well as numerically, only few things have been confirmed definitely in spite of huge 
number of studies. Thus, it is highly desirable to have analytic method for understanding 
the thermodynamic and critical behaviors in random spin systems. The method of gauge 
transformation is a useful tech niq ue to derive analytic results in gauge symmetric random 



spin systems (20 , 



21 



43. 



44. 



45 



46|,|47lQ- 

It was firstly applied to the Ising SG models 



and provides the internal energy and an upper bound on the specific heat exactly as non- 
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singular functions of the temperature on the Nishimori line in the temperature-randomness 
phase diagram. Further, it is proven that the boundary between the FM and non-FM 
phases in the temperature lower than the multicritical point (MCP, see Fig. at which 
paramagnetic (PM), FM and SG phases meet and the Nishimori line intersects, is vertical 
or re-entrant. Using a model with a modified probability distribution of randomness, the 
verticality is concluded Q|, which implies the absence of re-entrance: Although the theory 
needs a plausible but unproved assumption on the thermodynamic behavior of the modified 
model, the result is consistent with other theories and experiments of Ising-like spin glasses 



49]. The gauge transformation has been applied to gauge glasses with various symmetries 
20( including the U(l) symmetry, i.e. the XY gauge glass model. The same results as 
in the Ising case are derived. The absence of re-entrance is concluded for any dimensional 
model with a plausible assumption for the modified model, if the FM (or KT) phase is sta- 
ble in disordered regime. The gauge transformation is also extended to dynamical systems 
. For the Ising SG case, some important relations are derived for nonequi- 



mm 



47 



librium auto-correlation functions, order parameter and energy relaxations. They provide 
discussions on the absence of re-entrance, the regime of the mixed phase in the phase dia- 
gram, the multi-critical dynamics, and the equivalence of two nonequilibrium processes in 
the SG phase. 

The dynamics is one of most important aspects in disordered systems because of the 



slow relaxation In the case o: 

investi gate d for the mean field model 

HI54II55I] , for the aging phenomenon 
is a typical realization of slow dynamics, especially for SG materials 
contrast to the standard SG systems, not many studies on the dynamical properties have 
been performed for gauge glass systems. It would be helpful for future investigations to 
present exact results including dynamical properties. 

In this paper, we study the gauge transformation for dynamical systems of the gauge 
glass model. In the static case, the gauge symmetric models can be treated in the same 



standard SG theory, dynamical properties have been 



jjL for the Griffiths phase in short-range systems 



,0H 







I59I 6ol | and so on. The aging phenomenon 



62, 



63, 



64, 63. In 



manner as the Ising SG 



ory for the Ising SG 



45 



47 



2C I . It is also straightforward to apply the dynamical gauge the- 



481 ]. in which the stochastic dynamics is treated, to the gauge 



glass model. Since the gauge glass system has a continuous spin space, one can consider 
a deterministic dynamics with equations of motion, which is also treated in the present 
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theory. Several exact relations are derived among dynamical quantities such as equilibrium 
and nonequilibrium auto-correlation functions, relaxation functions of order parameter and 
internal energy. Using them, we discuss physical properties including the aging, the multi- 
critical dynamics, dynamics in the SG, a possible mixed phase and Griffiths phase, and the 
absence of re-entrant transition in two or higher dimensions. 

The organization of the paper is as follows. In section m the gauge glass model associated 
with the plane rotator is introduced. The probability distribution is chosen so that the 
system has the gauge symmetry. The results of the static gauge theory is summarized. In 
section a stochastic dynamics is introduced using the master equation formalism, and 
some dynamical quantities are defined. The gauge transformation is introduced for this 
system. Some exact relations are derived for dynamical quantities. In section HVl physical 
properties are discussed by use of the exact relations obtained. In section El equations of 
motion are introduced for the same gauge symmetric system, and the gauge transformation 
is applied. It is shown that the same relations can be derived for several choices of dynamics. 
The last section is devoted to the summary. 

II. GAUGE SYMMETRIC GAUGE GLASS MODEL 

The Hamiltonian we consider is 



where fa, uiij G [— n, n). The set <fi = (fa, fa, ■ ■ ■ , (fix) represents a configuration of total N 
spins, and u) represents a configuration of total N B quenched random variables associated 
to each bond. The summation of (ij) is taken over all bonds; we make no restrictions on 
the type or the dimensions of lattice, whereas one may suppose a usual nearest-neighbor 
interaction on the d- dimensional hypercubic lattice. The variable ujij has an odd property; 
Uji = —oJij [75J. For a particular bond configuration, the thermal equilibrium distribution 
at temperature T = J/k^,K is given by 






(2.2) 



The partition function is defined by 




(2.3) 
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where the integration of d<f> denotes the multiple integrations for the N spin variables (pi. 
The probability distribution for randomness is taken as 

P ^ D) ~ 11 2*I (D) ~ Y{D) ' (2 ' 4) 

w) 

where F = (0, 0, • ■ • , 0) represents an all-aligned spin state, I n (D) is the modified Bessel 
function, and 

Y(D) = (2nI (D)) NB . (2.5) 

In this distribution, which is central to the present paper (see also Fig. El below), the variable 
D controls the strength of randomness; D = and oo correspond to the full random case 
and the non-random case, respectively. The random average is denoted by 

dw P(u};D) ■■■ . (2.6) 

The integration of du; denotes the multiple integrations for the Nb random variables cjy. 
In the following, for simplicity, we will sometimes omit the dependence on spin set 4>, bond 
set oj, inverse temperature K and/or randomness D from functions defined above, if they 
are trivial or unimportant. 

The gauge transformations for functions of <f> and u) are defined by 

: (pi — > (pi - ipi (2.7) 
Vrf, : Uij — > LOij +ipi-ipj, (2.8) 

where tp = (ipi,ip2, • • • ,^n) is an arbitrary state of iV spins. Each set of transformation 
forms a group. While variables (pi and uJij possibly take values out of the domain, [— 7r, 7r), 
in this transformation, we always identify them with a residue of 2n. The Hamiltonian (j2.1j) 
is invariant under the transformation U^V^: 

U^Hifru) = (2.9) 

The gauge transformation of the distribution (|2.4jl is given by 



V^P(uj;D) = '- = Peq (i/>;D,u>), (2.10) 
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where we note that the spin set in Ti. and the temperature are different from the usual ones. 
Another important property in the theory is the invariant integral for (f> and uj: 



&4> ■■■ = J d<£ ••• , (2.11) 
do? ••• = J du V+ ■■■ (2.12) 

Using eqs. (|2.9|) and (J2] fUl, one can show [2^ that the partition function is gauge invariant 
under V^, 

V^Z(K;u) = Z{K;u), (2.13) 
providing the invariance of thermal distribution, 

Peq (<f>;K,u) = p eq ((f>;K,u>). (2.14) 

Then, we obtain the gauge invariance of the averaged energy; 

V* (H) k = (H)k. (2.15) 

Note that we will use the terminology "gauge invariant" only for functions of u> invariant 
under V^. Since the lhs of eq. (|2.12j) is independent of if), it is convenient to derive another 
expression of random average instead of eq. ([2.6)1 : 

[ "' ]c = / d " d ^ (JpTfDj P^ D ^ V * • ( 2 ' 16 ) 
The random average of gauge invariant function Q(u>) can be expressed as 

= / dw (Iffm < 2 ' 17) 

In three or higher dimensions, the topology of the K-D phase diagram is expected as in 
Fig.Q where PM, FM and SG phases appear. A possible Griffiths phase is also indicated. In 
two dimensions, the KT phase appears instead of the FM phase. The SG phase disappears 
or becomes a possible KT-like SG phase. The dashed line K = D is called the Nishimori-line 
around which following properties have been found from the static gauge theory [20! : 

• The energy and the upper bound of the specific heat are expressed by analytic functions 
of temperature; 

E(D,D) = -J^\nY(D) = ~N B J^, (2.18) 
C(D,D) < k B D 2 ^\nY(D). (2.19) 
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FIG. 1: Typical phase diagram of gauge glass model in the D-K plane {D controls the randomness 
and K = J/k-gT). The dashed line K = D is the Nishimori line. The FM phase boundary in low 
temperature region is assumed to be vertical. A possible regime of the Griffiths phase is indicated. 

• The FM and SG order parameters coincide with each other; 



Thus, no SG phase appears on this line, since m = in the SG phase. 

• The line is likely to intersect the MCP of PM, FM and SG phases. 

• The phase boundary between the FM (or KT) and the SG (or PM) phases below this 
line must be vertical to the D-axis or re-entrant j^j]. A plausible argnument can be 
made for the verticality. 

It is remarked that the probability distribution function (J2.4|) is different from the usual 
Gaussian form, which has been studied in many cases. This is because that we need a 
gauge symmetric form in the theory, which restricts the distribution function relating to the 
Hamiltonian as in eq. (j2.4j) . It is shown that the same properties can also be proven for 
another gauge symmetric U(l) model, the Villain model 66] with Gaussian random gauge, 
in which the local Boltzmann factor is given by 



m(D, D) = q(D,D). 



(2.20) 



oo 




(2.21) 
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FIG. 2: The probability distribution function for one bond p(x;D) = exp(D cos x) / (2ttIq(D)) . 

We consider that physical properties are equivalent in these gauge glass models irrespective 
of the detail of distribution P. In fact, as seen in Fig. the shape of the function ()2.4|) is 
quite similar to the Gaussian. As will be seen later, the present theory can be applied to 
the Villain model in the case of the stochastic dynamics in section ITTT1 However, we proceed 
our present theory by the rotator model. It is somehow difficult to apply in the case of the 
deterministic dynamics, since the explicit form of the Hamiltonian, which is necessary in the 
equations of motion, is not unique. 



III. STOCHASTIC DYNAMICS 



First, we consider the stochastic dynamics following the results of Ising SG (45 
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A. Master Equation 

We consider a Markov process for a fixed bond configuration, in which the density of 
state evolves with the master equation \(y\ . 
d 



dt 



Pt{4>) 



d# {W{4>\4t) ptW) - W{4>'\ct>) Pt {<t>)} . 



(3.1) 



The transition probability W(<p\4>') is non-negative and satisfies the condition of the detailed 
balance, 

W&Mp^MK^) = W(<l>'\<l>)p^frK,w), (3.2) 



S 



which guarantees the stability of the equilibrium distribution p eq . The condition (|3.2|) is 
automatically satisfied if one uses the expression with a symmetric matrix w(4>\<p'): 

= (3 ' 3) 



The function w{4>\4>') depends on the detail of dynamics considered. As in |45j], I assume 

W(4\4t) = ™0<W, 0'] Pc q (0) 0{A ^^Ve q (0') 0(A[0 ' 0D (3-4) 

for the single-spin-flip Metropolis dynamics j^, and 

for the single-spin-flip heat-bath dynamics, where Q(x) denotes the Heaviside function, 

A[0, ^=W(0;u;) w) (3.6) 
is the energy difference of two states, 

5 1 [^<f>'} = #1 (3.7) 

i 

^[0,01 = n ( 3 - s ) 

are the single-spin-flip and the zth-spin-flip operators, respectively. Note that the denomina- 
tor of the summand in eq. (j3.5)) is unchanged if the variable <p' in the argument of 5u [(f)", (ft'} 
is replaced by 4> providing the symmetry of w((f>\(p'). 

Since the diagonal element plays no role both in eqs. (j3.1|) and (|3.2jh one may define it 
in such a way as 

Then, the master equation ()3.1|) is reduced to 

lPt {cj>) = J tyt w{4>\4t) PtW), (3.io) 

and the solution is given by 

pt(4>) = J W <0|e* w |0') Po (0'). (3.H) 



dt' 
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FIG. 3: Illustrations of two nonequilibrium processes, I and II. Between t w and t + t w , the 
correlation is measured. 

The function (<j> |e tw | 0') is the integration kernel for time-evolution defined by 

00 j-n 

<0|e' w |0') = ]T-<</>|W"|0'>, (3.12) 

n=0 

// n— 1 \ n— 1 

ri d ^ (£) n w ^ (fc+i) i^ (fe) )' ( 3 - is ) 

\l=l / fc=0 

where (p^ = </> and 4>^ = <j>' . It has been proven that all eigenvalues of the operator W 
with (|3.2j) are real and negative except the zero eigenvalue corresponding to the equilibrium 
distribution j67|. Thus, all solutions with any initial conditions tend to the equilibrium 
distribution as t — > oo. 



B. Dynamical Processes 

We examine two nonequilibrium processes I and II in Fig. El The process I is related to 
the zero-field-cooling (ZFC). At the initial time t — 0, the system is kept in the equilibrium 
state with a temperature D and zero field; hereafter, we use a terminology "temperature" 
for K and D instead of "inverse temperature". The temperature is immediately changed 
(usually quenched) and the system relaxes in a heat bath with another temperature K in 



t > 0. The average for dynamical ensembles in t 



lis process is denoted by (•••)#■ The 



process II is related to the field quench (FQ) |59|,|63j]. The system starts from an all-aligned 
state F at t = and relaxes in the same heat bath as in the process I in t > 0; the average is 
denoted by (■ ■ • )J C . Since the all-aligned state provides the strong-field limit, this represents 
the process with a field quenched from oo to zero at t = 0. Note that the FQ is not equivalent 
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to the field cooling in which the applied field is weaker and is quenched after the waiting. 



We define the waiting-time dependence of nonequilibrium auto-correlation functions 



59 



601 ] for the processes I and II. In the language of master equation, that for the process I is 



expressed as 

^Ut+U)-iUU)^ D K = J d<Wd</> <0 |e <w | 0') (<f>' |e'- w | 0o> D) e^"^. 

(3.14) 

When K = D, the average {•••)% represents the dynamical average in the equilibrium 
process (• • • )^. In this case, the auto-correlation function ()3.14|) becomes independent of t w 
since the thermal distribution p e q{4>] K,lj) is an eigenfunction of W(4>'\4>), i.e. 

^ e i<j>i(t+U)-i<t>i(t^ D = ^gi^(t)-i^(0)^q (3.15) 

= J d<W (</>|e tw |<//) p^'-K)^-^. 

If one defines an auto-correlation function for the process relaxed from a fixed state ip - 
let us call this process the process II' — 

^ e i<Mt+tw)-i<MM^ = J d d 0' ^ | e *w| ^/^ | e t w W| ^ e i^-i^ ; (3 16 ) 

that for the process II is expressed as (V ( ^ ( * +iw - ) ~ 1( ^ ( '* w ' ) )^, and eq. 1)3.14)1 is rewritten as 

( e i^(t+*w)-i^(*w)^ = [dip p cq (^;D) (e^^J-^W)^ (3.17) 



In the long-time limit (t — > oo), the thermodynamic limit of ^e'^W-'^W^ ^j converges to 
the SG order parameter ^ea, which is non- vanishing in FM or SG phases 49, 50]. If one 
sets t w = 0, equation ()3.16j) with ip = F becomes 

(e^W>^ = J d0 <0 |e' w | F) (3.18) 

which is the FM order parameter at time t relaxed from the complete FM state, F, at 



approaches the spontaneous 



t = 0. When t — > oo, the thermodynamic limit of (e 1 ^*^)^. 
ordering. Note that all quantities defined above containing the exponential e 1 ^ have real 
values in the average because of the spin-flip symmetry, fa — > (pi + tt, for the Hamiltonian 
and the matrix W((f)\<p'). 
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Another relaxation functions, the nonequilibrium relaxations of exchange energy, are 
defined as 

{H{t)) D K = J j d<£d<£' <</>|e* w |0'> p eq (0'; A (3.19) 

(H(t))* K = jjdcj> (<f>\e tw \F) (3.20) 

C. Gauge Transformation 

Let us examine the gauge transformation of the functions defined above. We introduce 
the transformation for </>', 

% : (t>'i — ► ti-A, (3.21) 

to show the invariance 

U^U'^ {4>\e tw \<t>') = {<t>\e™\<t>'). (3.22) 

From eqs. ()2.9|) - (j2.13j) . it is easy to see that p e q(4>', K, u>), Su [<p, <fi'] and A[<p, <p'] are invariant 
under U^U^ V^. This yields the invariance 

U+U^Vi, w(<f>\<j>') = w{<p\<p') (3.23) 

in the case of the Metropolis (|3.4j) and the heat-bath dynamics (13 .5j) . It is natural to restrict 
the dynamics in such a way as eq. (|3.23|) holds. This yields the invariance, 

U^U'^W&W) = W{tp\<p'\ (3.24) 

where the invariant integral for 0" like eq. (j2.11|) is necessary to derive for the case (p — 4*' ■ 
Let us consider term by term in the expansion of ((f) |e* w | </>'), eq. (l3~T3j) . Similarly to 
and (f>', the gauge transformations for (f>^ (k = 1, 2, • • • , n — 1) are introduced as 

Uf : # } — 4$* ~ A- (3-25) 
Then, equation ()3.24|) with the invariant integral for <p^ like ()2.11|) leads to 

//n— 1 \ n— 1 

\i=l J k=0 

„ /n— 1 \ n— 1 

= y u w ^ (k+i) \^ ( 3 - 26 ) 



v^=l / k=0 
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providing the invariance (|3.22|) . 

From eqs. (|2.14|) and (|3.22|) with the invariant integral (|2.11j) taken into account, one can 
show that the auto-correlation function for the process I is gauge invariant as 



V. h (e 



i<fc(t)-i<fc(0)\eq 



)k = ( 
)l = < 



^i(t)-i^(0)\eq 



A'' 



1 i0i(t+tw)-i^i(tw)\- D 
/A-' 



(3.27) 
(3.28) 



while that for the process II is transformed as 



j4>i (t+t w )-i</>;(tw) N 



>>* = j d<W Ui,U'+V+ <</> |e tw | <//) (tf>' |e*- w | F> e^ 

_ ^ e i0i(t+tw)-i0i(tw)^ (3.29) 

where we use the invariant integral (j2.11|) . and 

^<0|e tw |0'> = ^<0|e' w |0') = <0|e tw |0' + V> (3.30) 

obtained from eq. (|3.22j) . The variables —ij) and (j)' + i[> represent the spin states, 
(—0i, —02, ■ • ■ , -tPn) and + Vi, 0' 2 + ^2, • • • , <#v + ^iv), respectively. Using eqs. (12. 1(1 j). 
(EHTjl and (|3~23jl . we derive 



(2n) N Y(D) V 



Because of the gauge invariance (|3.28|) . equation (|3.31|) is rewritten by eq. (|2.17|) as 



ji0i(t+t w )-i0i(<w) 



K 



(3.31) 



(3.32) 



for any waiting time £ w , any time interval £, any temperature K and any degree of random- 
ness D. We call it the "aging relation", since it relates aging phenomena in two distinct 
processes whatever waiting-time dependence is essential for the aging. From the invariance 
(12. 9|) . similar relation can be derived for the energy defined in eqs. ()3.19|) and ()3.20|) : 



(n{t)) 



(3.33) 



IV. PHYSICAL PROPERTIES 

The aging relation (j3.32j) and eq. ()3.33j) include parameters K and D characterizing the 
relaxation in t > 0. For each phase, one can examine the physical meaning of the relation by 
choosing (K, D) appropriately. Note that the initial temperature of the process I is always 
located on K = D. 
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A. Dynamics in the SG Phase 



First, let us consider the region with sufficiently strong randomness (small D) and suf- 
ficiently low temperatures (large K), where the SG phase appears if it exists. In such a 
region, D indicates a high temperature (small D), and the relaxations 



/ gi^i (i+iw ) -ifa (*w ) \ 



and 



H{t)) R for the process I describe a ZFC process (cooled immediately from a high 



3 i0i(t+t w )— i0i(£w) 



IK 



and 



for the 



temperature with D). On the other hand, 
process II are regarded as the relaxations in a FQ process (demagnetization from the strong 
field limit). Thus, relations ()3.32j) and ()3.33j) imply the equivalence of the dynamics es- 
pecially the aging phenomenon in the ZFC and the FQ processes. The same relation was 
pointed out in real SG materials (g^Hg^. At a glance, this relation expresses a trivial fact, 
since both initial states are located in the PM phase and both systems relax into the SG 
phase. However, the equivalence at any waiting time, which means the equivalence of the 
dynamical structure at any stage of relaxation, is non-trivial. 

In the case of sufficiently weak randomness (large D), the initial states of both nonequilib- 
rium processes represent almost the same spin state (complete FM state). Then, equations 
()3.32|) and ()3.33|) show a proper phenomenon. 



B. Dynamics on K = D 



Next, we examine the dynamics on the line K = D, where the dynamical average (■••}# 
becomes the average in the equilibrium process (• ■ ■ )^. From eqs. ()3.15|) and (|3.32J) . we 
derive 



^ e i<£i(t+tw)-i<M*w)^ p 



If one sets t w = 0, equation ()4.1|) yields 



^iW-i&COKeq 



[< 



^i(t)-%(0)\eq 



a, 



(4.1) 



(4.2) 



which implies the equivalence of FM and SG orderings (if one takes the limit of t — > oo). 
Since, in the SG phase, the SG order remains finite whereas the FM order disappears, eq. 
(j4.2j) indicates the absence of the SG phase on the line K = D as in the analysis of the 
static orderings p^Hjl Q- 
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Equations 1)4.1)) and (|4.2jl also imply the equivalence of the equilibrium and nonequilib- 
rium relaxations. The internal energy of the initial state F is identical with that on K = D, 
which is equal to the equilibrium energy on K = D (12.18)) [2^ . It is easy to see that, the 
rhs of eq. (J3.33)) is independent of time t if K = D; 

h(D) 



(4.3) 



Since the equilibrium energy is an increasing function of the temperature, we expect that 

F 

{"H(t)) K is an increasing function of t above K = D and a decreasing one below it. 

J c 

As seen from eqs. 1)4.1)) and (|4.2p . the auto-correlation function in the FQ process is 
independent of the waiting time, which suggests the absence of aging on K = D. In 
the region below the MCP on K = D, where the spontaneous ordering appears, it is not 
clear that the FQ process is appropriate for the observation of the aging, since the initial 
all-aligned state is not so different from the final equilibrium state with broken symmetry. 
While the rotational symmetry is initially broken in the FQ process, the dynamical behavior 
is equivalent to that in the equilibrium (symmetric) process — see eq. (14.1)1 . Thus, we assume 
that the FQ process exhibits an aging feature even in such a symmetry-broken region. 



Mixed Phase and Griffiths Phase 



The aging has been considered an inherent 
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63 



64, 



65j and theoretical 156 



57 



58, 
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)roperty in the SG phase from experimental 
60( view points. Since the aging phenomenon 



is a typical observation of the complex phase space for slow dynamics, it could occur in 
other complex phases. In such a phase, the aging would also be inherent, which means that 
it occurs whole of the phase when it is observed in some parts of the phase. 

Following the results of the SK model and experimental observations [49] for SG phe- 
nomenon, one may consider a possibility of the mixed phase between the FM and SG phases 
even in the gauge glass systems. It is natural to consider that the aging occurs in the mixed 
phase, since the SG feature in the mixed phase provides typical slow dynamics which reveals 
the aging. If the aging is inherent in the mixed phase, this indicates that the line K = D 
does not enter the mixed phase. This restricts the topology of the phase diagram: it must 
appear below K = D even if it exists. 

In the Ising SG systems, it has also been pointed out that there exists a dynamically 
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singular phase called the Griffiths phase 69] between the critical temperature of the pure 



system and the phase boundary of low temperature phase (the FM or the SG) 



53, 



5J, |55j - see Fig. ^ The possibility of the Griffiths phase has not yet been discussed 
much in the gauge glass systems. Figure Q is proposed just following an analogy with the 
dilute ferromagnet 69]. Concerning to the Griffiths phase, two cases can be considered to 
understand the present result if the aging is an inherent property: (a) There is no aging 
whole of the Griffiths phase, if the line K = D intersects the Griffiths phase as in Fig. 
(b) There is no Griffiths phase at least around K = D. In the former case, even if slow 
dynamics is observed in the Griffiths phase, it is quite different from that in the SG phase 
|70|. The latter case allows the existence of the Griffiths phase below K = D. Another 
possibility is that (c) the aging is not inherent at least in the Griffiths phase. This means 
that the aging occurs in some parts of the phase. Although it is not obvious which is correct 
in the present framework, the above result restricts the existence of the aging phenomenon 
and the region of the Griffiths phase, and supports future investigations of them. 

D. Multicritical Dynamics 

Let us consider the critical relaxation at the MCP. At the FM critical point, in general, 
the FM order parameter decays in a power-law asymptotically as 

m(t) = [(e^W>J c ~ r A ™. (4.4) 

Then, the dynamic scaling hypothesis [zi| 

m(t, e, L) = L-P /u rh(L 1/u e, L~ z t) e = (K c - K) /K c (4.5) 

reveals the relation of exponents 

A m = (4.6) 

zv 

On the other hand, the equilibrium auto-correlation function decays in a power-law as 

q(t) ee [( e ^(*)-^(o))^] c ~ r A « (4.7) 

at both FM and SG critical points. At the FM critical point, the exponent X q is associated 
with the FM criticality; a similar dynamic scaling reveals 

A, = ^. (4.8) 

zv 
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Thus, X m and X q are twice different in the FM case. At the SG critical point, it is associated 
with the SG criticality; the dynamic scaling hypothesis for the SG ordering 

q{t, e, L) = L'^qiL^e, L~H), (4.9) 

where tilde exponents are those for SG ordering, yields 

\ = JL (4.10) 

zv 

The scaling field is not unique in the random case. One can define another scaling form 
instead of ()4.5|) with another scaling field such as (D c — D)/D c . At the MCP, the scaling 
field e can be chosen both as e = (D c — D)/D c and as e = (K c — K) / K c since they are same 
order around the MCP on K = D. Note that critical exponents f3 and v depend on the way 
to approach the critical point. Thus, they are defined along K = D in eqs. (14. 5j) and (|4.fij) 
for the MCP case. On the other hand, the ratio (3/v and the exponent z are independent 
of the scaling field but depend on the critical point itself. On K = D, we have derived the 
relation (|4.2jl . This gives 

L-Pl v m{L xtv e, L~ z t) = L-^ D q{L 1/0 e, L~H), (4.11) 

for any L, t and e (K or D) providing 

A m = X q , (4.12) 

z = z, (4.13) 

P/v = P/v. (4.14) 



Note that the line K = D does not enter the SG phase but just touch it at the MCP 
The dynamical exponent z and the ratio fifv are independent of the way to approach the 
criticality, and characterized by the criticality itself. Therefore, exponents z and j3/u are 
associated with the SG ordering at the MCP. Equations (|4.12j) - (|4.14|) relate exponents of 
FM ordering and SG ordering at the MCP. The relations for FM dynamic exponents, ()4.6|) 
and (|4.8|h would keep along the phase boundary from the pure case up to the MCP in Fig. 
^ and so does the relation for SG ones (|4.10|) from D = up to the MCP; it is not clear 
whether the values of exponents are universal or not. 
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E. Absence of Re-entrance 



Here, we discuss the asymptotic behavior of the auto-correlation function for the process 



I with t, 



0. 



j<j>i(t)-i<t>i(0) 



which satisfies 



/ e i(j>i(t) -Wi (0)^ D 



K 



(4.15) 



to show the absence of re-entrant transition, i.e. the verticality of the FM boundary. We 
present the following two propositions for it. Since the long-time limit of this function is 
expected to behave a kind of order parameter, (i) the asymptotic behavior should be prescribed 
by both phases where the initial point (D, D) and the final point (K, D) locate. In Appendix 
lAl we mention that the average of gauge invariant function can not detect the FM boundary. 
Because of the gauge invariance (|3.27|) . (ii) the asymptotic behavior is independent of whether 
(D, D) and (K, D) are FM or not. On this point of view, the behavior of (q 1 ^)- 1 ^^)^ 

J c 

can change only at which (D, D) or (K, D) intersects the PM boundary. On the other hand, 
it approaches an FM order parameter because of eq. (|4.15|) . implying that the asymptotic 
behavior changes depending on if (K, D) is FM or not. At a glance, this contradicts with 
(ii) on the FM boundary below the PM boundary. However, they are not in conflict with 
each other if the line K = D intersects the MCP and the FM phase boundary is vertical in 
the temperature below it. In this case, when (K, D) locates on the FM boundary, (D, D) 
locates on the PM boundary simultaneously (i.e. the MCP). 

In two dimensions, the same argument can be made if one considers the asymptotic decay 
of relaxation functions as an order parameter for the KT phase instead of the long time limit 



of them. The order parameter relaxation 



^e 1 ^*))^. decays exponentially in the PM phase 

J c 

while it decays algebraically in the KT phase. Thus, the above argument is recognized as 
the statement that the KT boundary is vertical below the MCP. 



V. DETERMINISTIC DYNAMICS 

In section 11111 we treat the stochastic dynamics obeying the master equation with con- 
tinuous time (jH.ljh The results can also be applied to the discrete time interval case, which 
is used in Monte Carlo simulations. Since the spin variable is continuous in the gauge 
glass model, one can define equations of motion. In this section, we discuss the case of the 
deterministic dynamics like in molecular dynamics simulations. 
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It is noted that the essential features for a physical quantity Q to derive relations like 
(13.32)1 and ()3.33|) are the gauge invariance 

V+ (Q)k = (Q)k (5-1) 
for the process I as eq. ([3.28)1 . and the transformation rule 

V+ (Q)l = (Q)t (5.2) 
for the process II as eq. ([3.29)1 . By use of eqs. ([2.16)1 and ([2.17)1 . they reveal the relation 

[<Q)k)c = [(Q)k]c (5-3) 
irrespective of the detail of dynamics. 

A. Microcanonical Dynamics 

Let us introduce a conjugate momentum Li for each spin variable (pi and consider the 
total Hamiltonian Htot = J7~Ltot(L, 4>,u>) including the kinetic energy; 

7itat(L,0,a>) = - ^2 cos (^ ~ $i + + 7^^2 L2 i- ( 5 - 4 ) 

(ij) i 

The thermal distribution for this system is given by 
exp (-KH tot (L, 0,w)J 

— — — r = /i cq (L;iT) p eq (^;ir,u>), (5.5) 

Z tot {K;uj) 

where the partition function of 7i tot is defined by 

Z tot {K;u>) = jdLd(j) exp (-KH^, <!>,(*>)). (5.6) 

The integration of dL denotes the multiple integrations for iV momentum variables. For the 
spin variable, the thermal distribution of 7i tot is identical with that in the static case ([2.2)1 . 
The function // eq (L; K) is the thermal distribution for momentum variables; 

/ K \ N/2 K 
^ q (L;K) = (-) exp(- y ^L 4 2 ). (5.7) 

i 

The equations of motion for fa and Li are given by the canonical equations of Hamilton; 

dL d 

^ (jei) 

° n tot = JL t , (5.8b) 



dt dL 
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where the summation of (j G i) takes for all sites coupled with the % site. 

d d 

Note that the conjugate momentum Lj, differential operators — — and — — are invariant 
under gauge transformations and V^. Thus, the canonical equations (|5.8|) are also 
invariant under the gauge transformation U^V^. This means that if (L(t), 4>{t)) represents 
the solution for a bond configuration uj with an initial condition (Lo, 0o), the solution for the 
bond configuration V^u> with the initial condition (L , O — ip) is given by (L(t), <fi(t) — ■»/>). 
This implies the transformation rules 

u 'i, v i> <t>i{t] L o, 0o, w) = <fii(t; L , O , w) - V'ij ( 5 - 9a ) 
t/J, Lj(t; L , 0o, w) = L , 0o, w), (5.9b) 

where £7^ is the transformation operated to 0o- It is remarked that once the transformation 
for 0o is operated, it affects 0(t) at any time, since the solution is continuous in time. Thus, 
the operation U!^ in eqs. (I5.9J) can be recognized as U'^ in the stochastic dynamics. 

Let us define dynamical averages like eqs. (j3.14J) - ()3.16|) . Since the dynamics is determin- 
istic, the dynamical average is performed by the average over the initial condition (Lo, 0o)- 
When the initial distribution is in equilibrium, the temperature is unchanged since the total 
energy is conserved. The equilibrium dynamics with temperature K, ( • ' ' )^ — ( ■ • ' )^, is 
defined as 

{'"Tk = J dL o d 0o /U C q(L ;^) p eq (0 o ;^,^) (5.10) 

Although the dynamics contains no dissipative factor in this formulation, we assume that 
the system shows a relaxation to an orbit which traces the thermal distribution when one 
prepares an appropriate nonequilibrium initial distribution p o (L o ,0o)- This is relating to 
the ergodic theorem. There are several ways to determine po(Lo,0o)- One simplest way is 
to fix the total energy to the equilibrium value at the final temperature. In the process I, 
(•••)#, the initial distribution for spin variables is the thermal one with temperature D. 
The total energy is identical with the equilibrium value at temperature K providing 

Pl (L ,<j> ;K, D,u) = P ^'dZ) 6 (*M L o> 0o, «) - (H tot ) K ) . (5.11) 

where (H tQ t) K = N/2K + (H) K , and 

Qi(K,D,ut) = J dL o d0 o Pcq (0o; A^) 5 (w t ot(L O) 0o,u>) - (H tot ) K ) (5.12) 
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is the volume of the phase space in the process I. When K = D, the dynamical average with p\ 
becomes identical with (j5.10|) . because of the equivalence of the canonical and microcanonical 
ensembles in large enough systems. Similarly, the initial distribution for the process II' is 
given by 

p^(L ,(j> ;K,u>) = y° ~ ^ 6 (n m (L , cf> , u>) - (H tot ) K ) . (5.13) 

where 

= JdL d<f> 5(<f> -^)5(HUlo,<f> ,u)-(Htot)K) (5.14) 

is the volume of the phase space in this process. Using these initial distributions, we define 
the correlation functions for processes I and II; 

( e i*(t+tw)-i*(tw))£ = J dL o d0 o pi(L ,(j>o;K,D,uj) ^Ht+^-M^^Mo^) ^ 

(5.15) 

^ e %(t+tw)-^(tw)^ = J dLod^o p F (L o ,0o;^,w) e*^^'*"")-*^^'*""). (5.16) 

The dynamical averages for the energy are also defined in a similar way. 

Noting the invariance ()2.14j) - (j2.15j) . the function pi(Lo, <f>o] K, D, u)) is invariant under 
the gauge transformation U'^V^,; 

Pi {L ,(I>o;K,D,uj) = Pl (L ,<t> ;K,D,u). (5.17) 

This reveals the gauge invariance process I as in the stochastic 

dynamics case ([3.28)1 . On the other hand, the transformation for p-p gives 

UipVyp p F (L ,(f) ;K,u>) = p^(L o ,0 o ;^,^) (5.18) 

which reveals the transformation rule for the process II, 

Vj, {e^i{t+u)-\4>i{u)\^ = J dL d</> p^(L ,4>o;K,uj) e i ^ t+ ^ ^°^- i ^ u '^ 

Thus, form eqs. (|5.1|) - ([5.3|h equation (|5.19|) yields eq. ()3.31|) providing the aging relation 
f!3.32|) . It is easy to see that eq. ()3.33p also holds. 
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B. Liouville Equation 



To compare the formulations in stochastic and deterministic dynamics discussed above, it 
is instructive to consider the classical Liouville equation which describes the time evolution 
of distribution function p t under the canonical equations ([5.8)1 : 

dp t _ ^ f d^tot dp t dUtot dp t \ . . 

dt \ dU d& dfa dU J ' 1 • j 

If we define the solution of eq. (|5.2U|) for the initial condition p = 5(L — h )S(cp — cp ) as 
G(t, L, 0; L , 0o, w ), this is nothing but a delta function 

G(t, L, cp; L , 0o, w) = 5 (L - L(t; L , <p , w)) 5 (0 - <£(t; L , O , w)) , (5.21) 

since the orbit is deterministic. The function G is called the Green function, and the solution 
for any initial distribution po(L, cp) is given by 

p t (L,(p;u}, [po]) = f dL d(po G(t,L, </>;L , 0o,«) Po(Lo,0o), (5.22) 

where p f is a functional of the initial distribution indicated in the arguments as [po]. 

The Green function G plays a similar role with (cp |e iW | 0') in the stochastic dynamics. 
It is easy to see that eqs. (|5.9|) and (j5.21|) yield the invariance, 

U'j, G(t, L, cp; L , cp ; u) = G(t, L, cp; L , O ; w), (5.23) 

which is the same relation as eq. ()3.22|) . Since the differential operators and 7i tot is invariant 
under U^V^, when the function p t is the solution of eq. (|5.20J) with an initial condition 
Po, the function p t is also the solution with the initial condition po : It is easily 

confirmed from eq. ()5.22j) that 

UrpVj, pt{L,cp; u, [po]) = p t {L,cp;uj,[U^p }). (5.24) 

Using the definitions above with the Green function, the auto-correlation functions are 
expressed as 



(eWt+twHfcCtw)^ = J dLd</>dL'd<//dL o d0o (5.25) 
G(t, L, cp; V, cp', w) G(t w , V, cp'; L , <p , w) pi(L , O ; £>, w) e 1 ^, 
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for the process I, and 

j<pi(t+t w )-i<j>i(t w )\ F 

G(t, L, 0; V, 0', u) G(t w , V, 4>'- L , 0o, w) p F (L , o ; w) e 1 ^"^. 



= I dLd0dL'd0'dL o d0 o (5.26) 



for the process II. They are equivalent with eqs. (|5.15|) and ()5.16|) . From invariance ()5.23|) 
and ()5.17|) with eq. (j2.11j) . one finds that ^e 1< ^ t+ * w - )_1( ^ tw ' ) )^ is gauge invariant and the other 
one obeys 

^ ( e i*(t+M-i*(*w)^ = ydLd^dL'd^'dLod^o (5.27) 
G{t, L, 0; L', cj>', ui) G(t w , V, 0'; L , 0o, w) p^(L , O ; «) e^"^. 
Thus, form eqs. (|5.1|) - (j5.3|) . this yields eq. (|3.31|) providing the aging relation (|3.32|) . 

C. Canonical Initial Distribution 

Another way to define the initial distributions instead of (j5.11|) and (j5.13|) is given by 
the canonical ensemble. In both processes, the initial distribution is determined so that the 
total energy distributes in terms of the thermal distribution with temperature K; 

/T> ( E- u) e~ KE / J 
dL d0 o Po(L , 0o) S (E - ^ tot (L , O , w)) = \ I , ■ (5.28) 

where V tot (E; u>) is the density of state; 

V tot (E; w) = J dL d0 o <f(£ - H tot ). (5.29) 

For the process I, po has a form of 

Pi(Lq, 0o) = Pi(Lo,0o;^,A^) Pcq(0o;A^)- (5.30) 

Equation (|5.28j) implicitly defines the function pi(Lo, 0q; K, D, u>) depending on K, D and 
u>. For the case K = D, pi = p eq (L ; D) satisfies eq. ()5.28|) . For the process IF, po has a 
form of 

PV>( L o,0o) = Pv( L o,0o;^,^) ^(00 - tp)- (5.31) 

Equation ()5.28|) implicitly defines the function p^,(L , 0o; K, u>) depending on K and u>. 

Noting the invariance ()2.9|) . ()2.13|) and ()2.14|) with the invariant integral ()2.11|) . one finds 
that, if po(Lo, 0o) satisfies eq. ()5.28|) . so does ULV^, po(Lo, 0o)- Since U'^V^ pi instead of pi 
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satisfies the same equation as ()5.28j) with (|5.30p . the function pi(L , (po) is invariant under 
U^Vip. This reveals the gauge invariance of ^e 1 ^ i( -* + * w - )_1 ^ i( ' iw - ) )^ for the process I. On the 
other hand, if Pf(L , <p ) satisfies eq. (|5.28|) for the process II, so does 

p F (L , (f> ) = [U 1 ^ p F (L , cp ; K, u)] 5(</>o - V0- (5.32) 

This means 

U'^Vj, p F (L ,<t>o;K,u>) = p^(L ,(f} Q ;K,uj). (5.33) 

Therefore one can derive the transformation rules like (j5.19|) providing ()3.3H) and the aging 
relation (J3.32)) for the canonical case. 



D. Extended System 

In the above, we have defined relaxation processes in conservative dynamics with gauge 
symmetry. Here we consider another deterministic dynamics, which keeps the temperature 

d the extended system dynamics. 



constant 



In addition to the momentum variables, we introduce a single degree of freedom, the 
daemon, which couples to the original system. It will act as a heat bath for the original 
system and convert the dynamics as an fixed-temperature one. Let the variable y be the 
daemon and i the conjugate momentum of it. The Hamiltonian of the extended system 
H cxt = J7^ cxt (L, 4>, w) is defined by 



H cxt (L,cp,u>) = - ^ cos(& - 4>j + Uij) + Yl L i + ~2 + k V > (5 ' 34) 

(a) y i 

where <f> represents the set of spin variables (p including the daemon y, and L represents 
the set of momentum variables L including £. Note that the temperature K is explicitly 
included in Ti ex t- The canonical equations are given by 

= - J Y sin (^ ~~ fa + w »j")> (5.35a) 
^ = JUlv\ (5.35b) 

% = J ( £ L "/y 2 - ( N + iv> ( 5 - 35c ) 

^ = JL (5.35d) 
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It is known that the microcanonical ensemble of TC ex t with the total energy E becomes the 
canonical ensemble of 7i to t with temperature K irrespective of E when one traces out the 
daemon variables y and I. This means that the solution of the equations with an arbitrary 
initial condition shows a relaxation to an orbit which traces the equilibrium distribution 
with temperature K, if one assumes the ergodic theorem. 
In this system, the dynamical averages are given by 

(■■■)k = J d £od0 o fJ.(L ,£,y) p eq (0 o ;A^) (5.36) 
for the process I and 

(•••)* = JdL d0 o /i(L , *, V) <K4>o - F) • • • , (5.37) 

for the process II, where the integrations of dLo and d0o denote the multiple integrations 
for N + 1 momentum variables including the daemon's i and those for iV + 1 spin variables 
including y, respectively. As stated above, the initial distribution function /i(Lo, £, y) for 
L, £ and y can be chosen arbitrary. We assume that /i(L ,£, y) is independent of uj so 
that it is gauge invariant. It is natural to consider that the daemon variables y and i are 
invariant under the gauge transformation U^V^p. Thus, the canonical equations ()5.35|) are 
also invariant under the transformation, which reveals the transformation rule like (|5.9j) : 

U^p V$ (j)i(t; L , 0o, w) = <j)i(t; L , O , w) - ipi, (5.38) 

Following the previous subsection, one can show the gauge invariance of ^e 1 ^( t+t ^)- 1 ^( t ^)^ 
for the process I and the transformation rule like 1)5.19)1 for the process II providing (j3.31j) 
and the aging-relation ()3.32|1 for the extended system. 



VI. SUMMARY 



We have applied the gauge transformation to dynamical systems of the gauge glass. The 
stochastic dynamics with the master equation formalism (|3.1j) as well as the deterministic 
one with the canonical equations of Hamilton are treated. For the deterministic case, we 
consider standard equations of motion with no dissipative factor ()5.8)1 . in which one realizes 
relaxation processes to choose appropriate initial distributions. Both microcanonical and 
canonical ensembles are considered for the initial distribution. Further, we consider the 
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extended system dynamics ()5.35|) in which the temperature is fixed. In this dynamics, a 
single degree of freedom is introduced, which is called the daemon and plays a role of a heat 
bath. For all dynamics, we show the theory of gauge transformation in a coherent manner. 

The dynamical relations (|3.32j) and (|3.33|) are derived exactly for all dynamics. The 
equivalence of nonequilibrium relaxations in the SG phase has been shown between the ZFC 
and FQ processes. The waiting-time dependence, which is typical in the aging-phenomenon 
coincides for ZFC and FQ processes. On the line K = D,we have found that the equilibrium 
relaxation coincides with the nonequilibrium one from the strong field limit. This indicates 
the absence of aging providing the restrictions for the regime of possible mixed phase and 
Griffiths phase. The exact relations for critical exponents are found at the MCP on this 
line. Further, we have confirmed the verticality of the FM phase boundary or equivalently 
the absence of the re-entrant transition. 

In contrast to the standard SG systems, not many studies on the dynamical properties 
have been performed for gauge glass systems. The present results would be helpful for future 
investigations including dynamics. 
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APPENDIX A: MODIFIED MODEL 



The modified model 2JJ, |44| is useful to consider the qualitative behaviors of physical 
quantities and the topology of phase diagram. The Hamiltonian is identical with the original 
model. The distribution of randomness, including a fixed constant a, is modified from the 
original model (a = case), 

Y(D + a) Z(D;uj) 



P(u:D + a) 



Y(D) Z(D + a;u) 
The random average in the modified model with a is denoted by 



(Al) 



(A2) 
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It is shown |2C| that the average of any gauge invariant quantity is independent of a, i.e. 



{Q(u>)} a c = [Q(u,)] t 



(A3) 



In this appendix, we summarize the properties of the modified model obtained previously. 
Analyzing the FM and SG order parameters defined from the static FM and SG correlation 
functions, 



(A4) 
(A5) 



m a (K } Df = lim {(e^)^ 
q a {K,Df = lim {\{^~^) K 

the following properties have been found exactly for the model with a: 

• The PM boundary, at which the edge of q a = locates, is unchanged with a. 

• The ordered phase on K = D + a must be the FM phase, since m a = q a holds on it. 

• The line K = D + a is likely to intersects the MCP. 

• The FM boundary below K = D + a is vertical or re-entrant. Further, the verticality 
can be shown if the ordered phase between the line K = D + a and the non-random 
case (D = oo) is always FM for the modified model with any a; this assumption is 
quite plausible since both sides exhibit only PM-FM transition. 

These properties are identical with those for the original model (a = 0). In two dimensions, 
one can treat KT transitions in the same way by analyzing the correlation length for FM 
and SG orderings, 

-l 



ia(K,D) 



lim 



lim 



d 



dRij 
d 



In 



dlL 



lnjKe^-^) 



K I 



(A6) 
(A7) 



instead of order parameters (jA4j) and (jA5|) . 

Since the FM correlation function as well as the FM and KT order parameters eqs. (|A4|) 
and (jA6|) are not gauge invariant, the FM (KT) boundary changes with a: At least on 
K = D + a, the FM (KT) regime changes with a. Thus, the qualitative behavior of averaged 
gauge invariant quantities at (K, D) should not be influenced by the fact whether (K, D) 
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is FM (KT) or not. In other words, any gauge invariant quantity can not be an order 
parameter for the FM (KT) phase in modified models including the original one (a = 0). 
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